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ABSTRACT 

Heavily obscured active galactic nuclei (AGNs) play an important role in con- 
tributing to the cosmic X-ray background (CXRB). However, the AGNs found in 
deep X-ray surveys are often too weak to allow direct measurement of the column 
density of obscuring matter. One method adopted in recent years to identify heavily 
obscured, Compton-thick AGNs under such circumstances is to use the observed mid- 
infrared to X-ray luminosity ratio as a proxy for the column density. This is based on 
the supposition that the amount of energy lost by the illuminating X-ray continuum to 
the obscuring matter and reprocessed into infrared emission is directly related to the 
column density and that the proxy is not sensitive to other physical parameters of the 
system (aside from contamination by dust emission from, for example, star-forming 
regions). Using Monte Carlo simulations, we find that the energy losses experienced 
by the illuminating X-ray continuum in the obscuring matter are far more sensitive to 
the shape of the X-ray continuum and to the covering factor of the X-ray reprocessor 
than they are to the column density of the material. Specifically we find that it is 
possible for the infrared to X-ray luminosity ratio for a Compton-thin source to be 
just as large as that for a Compton-thick source even without any contamination from 
dust. Since the intrinsic X-ray continuum and covering factor of the reprocessor are 
poorly constrained from deep X-ray survey data, we conclude that the mid-infrared 
to X-ray luminosity ratio is not a reliable proxy for the column density of obscuring 
matter in AGNs even when there is no other contribution to the mid-infrared lumi- 
nosity aside from X-ray reprocessing. This conclusion is independent of the geometry 
of the obscuring matter. 
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1 INTRODUCTION 

Although it has been known since the 1980s that a significant fraction of the active galactic nuclei (AGNs) that contribute to 
the CXRB X-ray background (CXRB) must be heavily absorbed (Setti & Woltjer 1989), the details of the absorbed population 
remain uncertain (e.g., Comastri et al. 1995; Tozzi et al. 2006; Polletta et al. 2006; Gilli, Comastri, & Hasinger 2007; Frontera 
et al. 2007; Martfnez-Sansigre et al. 2007; Ballantyne & Papovich 2009; Treister, Urry, & Virani 2009; Draper & Ballantyne 
2009; Fiore et al. 2009; Georgantopoulos et al. 2008, 2009; Malizia et al. 2009, 2010). One of the reasons for this is that the 
problem is highly degenerate. Another reason is that, although much progress has been made in the last decade or so in 
resolving the hard X-ray background below 10 keV, the bulk of the resolved AGNs are too weak to robustly measure their 
column densities with X-ray spectroscopy. As a result, a number of indirect "proxies" have come into use for interpreting 
deep survey data and identifying heavily obscured AGNs. One of these is the use of the equivalent width (EW) of the Fe Ka 
line at ~ 6.4 keV as a rough indicator of the column density because EWs larger than ~ 1 keV cannot be produced by 
a Compton-thin X-ray reprocessor (e.g., see Ghisellini, Haardt, & Matt 1994; Ikeda, Awaki, & Terashima 2009; Murphy & 
Yaqoob 2009, hereafter MY09). However, AGNs found in deep surveys are seldom bright enough to allow a direct robust 
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measurement of the Fe Ka line EW. An additional constraint could in principle be obtained from an independent indicator 
of the intrinsic AGN continuum luminosity such as an emission line in a different waveband. A few have been explored, the 
ratio of the luminosity of [O in] A 5007A luminosity to the observed X-ray luminosity being the oldest (e.g., Maiolino et al. 
1998; Bassani et al. 1999; Zakamska et al. 2003; Cappi et al. 2006; Panessa et al. 2006; Reyes et al. 2008; Bongiorno et al. 
2010). Melendez et al. (2008) have explored the use of the ratio of the luminosity of the [O vi] A 25.89/im emission line to the 
observed X-ray luminosity as a diagnostic of the column density. Gilli et al. (2010) have recently suggested the analogous use of 
the [Ne v] A 3426A emission line. However, derivation of the column density from these indicators is highly model-dependent 
as there are degenerate solutions (in addition to uncertainties in calibration of the relations), and even then, the method 
can only give a line-of-sight column density, which may be different to that out of the line-of-sight. A third indicator that 
is used to identify heavily obscured AGNs in deep X-ray surveys is the ratio of the mid-infrared luminosity to the observed 
X-ray luminosity (e.g., Martmez-Sansigre et al. 2005; Daddi et al. 2007; Hickox et al. 2007; Alexander et al. 2008; Fiore et 
al. 2008, 2009; Georgantopoulos et al. 2008, 2009, 2010; Vignali et al. 2010; Bauer et al. 2010; Georgakakis et al. 2010, and 
references therein). The idea is that the intrinsic AGN continuum that intercepts the obscuring matter loses energy due to 
the absorption and Compton-scattering of X-ray photons and that this energy loss ultimately heats the medium so that it 
is then reradiated in the infrared band. It is then supposed that the ratio of the infrared to observed X-ray luminosity in 
AGNs could be an indicator of the column density of the reprocessor. Since reprocessed emission from dust can also provide 
a significant contribution to the infrared band (e.g., Maiolino et al. 2007; Horst et al. 2008; Rowan-Robinson, Valtchanov, 
& Nandra 2009; Georgakakis et al. 2010), it has certainly been recognized that interpretation of the mid- infrared to X-ray 
luminosity ratio indicator is complex. Nevertheless, it is still believed (e.g., Mullaney et al. 2010, and references therein) that 
it will eventually be possible to discriminate between different scenarios and utilize the mid-infrared to X-ray luminosity ratio 
to identify Compton-thick AGNs. In the meantime, whilst an unusually large mid-infrared to X-ray luminosity ratio does not 
necessarily mean that an AGN is heavily obscured, an unusually small ratio could (it is supposed) be used to establish that 
an AGN is not a candidate to be a Compton-thick source. In other words, rejection of Compton-thick candidates may also be 
useful. 

Unfortunately it has never been demonstrated, even without the complications arising from dust emission, that the mid- 
infrared emission to X-ray luminosity ratio should be sensitive to the X-ray absorbing column density to the extent that the 
ratio could used to infer the magnitude of that column density. One cannot correct for the effects of dust if the relation between 
the X-ray energy reprocessed into the infrared band has not been quantified theoretically. It is the aim of the present paper to 
investigate and quantify the relationship between the energy lost by the intrinsic X-ray continuum in the obscuring matter and 
the column density of that reprocessor, explicitly without the complications of dust emission. Calculating the energy available 
from X-ray losses for reemission in the mid-infrared band is of course only the first step in the process. However, establishing 
that the energy budget from X-ray losses is sensitive to the column density, and not sensitive to other key parameters, is a 
minimal requirement for the mid-infrared to X-ray luminosity ratio to be a useful indicator of the column density. If we find, 
as we do, that this minimal requirement is not met, then additional physical processes can only lessen the usefulness of the 
diagnostic. Our investigation is based upon the toroidal X-ray reprocessor model of MY09, but our principal conclusions will 
have general applicability and will be independent of geometry. 

The paper is organized as follows. In §2 we give a brief overview of the assumptions of the MY09 model and describe 
the method we used to analyze the energy budget. In §3 we present the results of X-ray energy losses in the reprocessor as 
a function of injection energy, and in §4 we present the results of integrating these "loss" spectra over an incident power-law 
continuum. It is in §4 that we show the relationship between the integrated energy loss and the column density of the obscuring 
matter, as well as the dependence of the energy losses on the covering factor and on the power-law spectral index. In §5 we 
define a proxy for the ratio of the infrared luminosity due to X-ray energy losses, to observed X-ray luminosity and present the 
results of computations of this ratio as a function of the column density of obscuring matter. We summarize our conclusions 
in §6. 



2 METHOD 

Here we give a brief overview of the critical assumptions that our model Monte Carlo simulations are based upon. Full 
details can be found in MY09. Our geometry is an azimuthally-symmetric doughnut-like torus with a circular cross-section, 
characterized by only two parameters, namely the half-opening angle, do, and the equatorial column density, Nu (see Fig. 1 in 
MY09). We assume that the X-ray source is located at the center of the torus and emits isotropically and that the reprocessing 
material is uniform and essentially neutral (cold). For illumination by an X-ray source that is emitting isotropically, the mean 
column density, integrated over all incident angles of rays through the torus, is Nh = (n/4:)Nu. The inclination angle between 
the observer's line of sight and the symmetry axis of the torus is given by o h s , where 9 b B = 0° corresponds to a face-on 
observing angle and 9 b B = 90° corresponds to an edge-on observing angle. In our calculations we distribute the emergent 
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photons in 10 angle bins between 0° and 90° that have equal widths in cos# obs , and reler to the lace-on bin as #1, and the 
edge-on bin as #10 (see Table 1 in MY09). 

We used a version ol our Monte Carlo code that injects single photons (as opposed to bundles ol photons, with weights), 
so that energy is explicitly conserved in the code. Photons were injected with energies in the range 5-500 keV. In the present 
study we are interested in the energy losses and not the emergent spectra so lor energies below 5 keV we calculated energy 
losses lor absorption only, since the relative contribution from Compton scattering is negligible in this regime. II the solid angle 
subtended by the torus at the X-ray source is AO, the covering lactor is [Afi/(47r)] = 1 — cos6* bs (#obs only varies between 
0° and 90° ) . The covering lactor may also be expressed in terms ol the physical dimensions ol the torus. II a is the radius 
ol the circular cross-section ol the torus, and c + a is the equatorial radius ol the torus then [Af2/(47r)] = (a/c) (see MY09). 
Our model employs a lull relativistic treatment ol Compton scattering, using the lull differential and total Klein-Nishina 
Compton-scattering cross-sections. We utilized photoelectric absorption cross-sections lor 30 elements as described in Verner 
& Yakovlev (1995) and Verner et al. (1996) and we used Anders and Grevesse (1989) elemental cosmic abundances in our 
calculations. The Fe Ka, Fe K/3, and Ni Ka fluorescent emission lines were included in the model, as described in MY09. We 
performed Monte Carlo simulations for (a/c) in the range 0.1-1.0 and Na in the range 3 x 10 22 cm~ 2 to 10 25 cm~ 2 . 



3 ENERGY LOSS SPECTRA 

For photons injected into the torus with an energy E, in the interval dE, we define D[(a/c), Nh, E]dE as the fractional 
energy lost by those photons due to X-ray absorption and Compton downscattering in the torus. This is shown symbolically 
in equation 1, where F a bsorbcd(-EO and Fc(E) are fractions ol the original energy E that are ultimately absorbed or lost to 
the medium through Compton downscattering respectively. iVo is the number ol photons injected but obviously cancels out 
ol the definition. The frmction D does not ol course only depend on (a/c), Nh, and E, it depends on many other lactors such 
as element abundances and geometry. A photon that undergoes Compton scattering is ultimately either absorbed or escapes 
the medium. However, absorption can result in line emission and such line photons either escape (Fl{E) in equation 1) or are 
themselves absorbed. 



D[(a/c),N H ,E] = N E[F^^(E) + Fc(E)-F L (E)] 
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In practice, using our Monte Carlo code we can track photons injected in a given energy interval and count the energy 
that escapes. The energy that does not escape must have been given up to the reprocessing medium and this is represented 
in the second part ol equation 1. Whereas the lunctions -F a b sor bcd(-E), Fc(E), and Fl{E) are already integrated over the 
energy distribution resulting from a single injection energy, E, the frmction F CBCa , pc (E, E') in equation 1 needs to be explicitly 
integrated over the distribution ol downscattered energies, from to E. We also summed over all escaping angles since we 
are interested in the total energy lost to the medium. Our model includes only three emission lines so we miss some ol the 
escaping energy and therefore overestimate the energy loss, which is already just an upper limit on the total energy available 
for reprocessing into infrared emission. Moreover, our conclusions in the present paper will not be affected by neglect ol 
the energy in the omitted emission lines. Note that il the torus were completely opaque and non-reflective at all energies, 
D[(a/c), Nh, E] would simply be equal to (a/c), the covering lactor, since that is the fraction ol the energy that would be 
intercepted and lost to the torus. 

In Fig. 1 we show curves olthe quantity D[(a/c), Nh, E] /(a/c) for five different column densities in range 3x 10 22 cm -2 to 
10 25 cm -2 (solid, black curves), computed from our Monte Carlo simulations. The curves have been calculated for (a/c) = 0.5 
but have then been normalized by the covering lactor, (a/c). The purpose ol this is to later lacilitate direct comparison ol the 
shapes ol the curves for different values ol (a/c). The units for the curves are fractional energy loss per unit energy, per unit 
covering lactor. We see that at low energies, below ~ 1 keV, the fractional energy loss per keV has already reached its limiting 
value ol (a/c) for all ol the column densities shown. At higher injection energies less ol the energy is lost but the threshold lor 
it to drop below (a/c) depends on Nh. For the two highest values ol Nn shown (5 x 10 24 cm~ 2 and 10 25 cm -2 ), this energy 
loss drops below the maximum only above ~ 10 keV. The drop in the curves in Fig. 1 due to the decrease in absorption with 
increasing photon injection energy for each JVh is followed by a flattening again and then a rise as the energy losses due to 
Compton downscattering become important. For the two lowest column densities in Fig. 1 we have illustrated how Compton 
downscattering losses compare with absorption by overlaying the equivalent curves for pure absorption only (red curves). 

We found that the normalized lunctions D[(a/c), Nh , E]/(a/c) could not be distinguished in magnitude or shape (within 
the statistical errors ol the Monte Carlo results) for any value (a/c) in the range 0.1-1.0 lor column densities up to 10 24 cm -2 . 
For higher column densities a dependence ol the magnitude and shape ol the fractional energy loss curves on the covering 
lactor does become apparent. This is shown in Fig. 2, which shows the curves in Fig. 1 lor Nh — 5 x 10 24 cm~ 2 and 10 25 cm~ 2 , 
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Figure 1. Curves of D[(a/c), Nh , E] (see equation 1), showing the difference between the energy escaping the torus and incident energy 
(per keV), as a fraction of the incident energy (per keV). The horizontal axis corresponds to the incident energy. The curves were 
calculated for a covering factor of [Afi/(47r)] = (a/c) = 0.5 but they are normalized by the covering factor so that they can be directly 
compared with curves with different covering factors (e.g., see Fig. 2). The black curves show the Monte Carlo results for five different 
column densities, Nh, as shown. The red curves show, for comparison, the effect of neglecting Compton scattering (i.e. for the case of 
absorption only) for the two lowest column densities. 

zoomed in on the region in which Compton scattering dominates. The black curves show the previous results with (a/c) = 0.5 
and the red and brown curves correspond to (a/c) = 0.1 and (a/c) = 1.0 respectively. The dependence on covering factor in 
the Compton-thick regime arises because the amount of energy reflected from the inner surface of the torus and intercepted 
by the torus again depends on the opening angle of the torus and therefore on the covering factor. 

The fractional energy loss curves in Fig. 1 and Fig. 2 do not depend on the shape of the incident X-ray continuum 
spectrum because the curves give the fractional energy loss for monoenergetic photon injection. In order to calculate the 
absolute energy losses the curves must be multiplied by the incident continuum energy spectrum, as in equation 2 (where A is 
an arbitrary flux normalization). Since the intrinsic X-ray continuum of AGNs can typically be characterized by a power law 
with a photon index, T, in the range 1.5-2.5, we can expect that the energy losses due to Compton downscattering compared 
to absorption will be diminished more for steeper spectra than for flatter spectra. 



L[(a/c),N H ,E] = D[(a/c),N H ,E]AE~ 



keV cm 



kcV~ 



(2) 



4 CUMULATIVE AND TOTAL ENERGY LOST 

So far we have examined the energy losses for monoenergetic injection. We now examine the energy losses integrated over a 
range of energy. We define the cumulative energy loss, CL[(a/c), Nh ,T, E] for a power-law incident photon spectrum of the 
form AE~ T , between a lower energy, El and an energy E, as a fraction of the total incident energy in the energy range El 
to some upper energy, Eu- Thus, we have 



CL[(a/c),N H ,r,E] = 



f* L L[(a/c),N H ,E} dE 
\** u iFr+i dE 



(3) 



The absolute normalization of the incident continuum of course cancels out in equation 3. In the remainder of this paper 
we use El = 0.5 keV and Eu = 500 keV. The cumulative fractional energy loss function as defined in equation 3 is useful 
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Figure 2. The dependence of the energy loss functions, D[(a/c), Nh, E], on the covering factor, AS) / (Air) = (a/c). Shown are curves 
for the two highest column densities in Fig. 1, Ng = 5 X 10 24 cm -2 and 10 25 cm" 2 . The dependence on (a/c) for Njj = 10 24 an" 2 is 
negligible compared to the statistical errors. Each curve has been normalized to the particular value of (a/c) for that curve in order to 
facilitate comparison of the shape of the functions. The values of the covering factor, (a/c), in both (a) and (b) are 0.1 (red), 0.5 (black), 
and 1.0 (brown). 



for showing the energy range over which energy losses are most important (for a given set of the parameters (a/c), Nh, and 
r). In Fig. 3 we show calculations of the function CL[(a/c), Nh, T, E] in equation 3 for a covering factor of 0.5 and three 
different values of F (1.5, 1.9, and 2.5), each for two different values of Nh (red curves correspond to 3 x 10 22 cm -2 and black 
curves correspond to 10 25 cm -2 ). These two extreme column densities correspond to equatorial Thomson depths of the torus 
of ~ 0.024 and ~ 8.1 for the smaller and larger column respectively, and are therefore representative of a Compton-thin and 
a Compton-thick case respectively. In Fig. 3, the total energy loss in the range 0.5-500 keV as a fraction of the total incident 
energy in the same energy band can be read off from the curves at E = 500 keV. 

Two results are immediately apparent from Fig. 3. The first is that changing the column density by a factor of more than 
300, from 3 x 10 22 cm -2 to 10 25 cm -2 produces a disproportionately small change in the total integrated fractional energy 
loss, CL[(a/c), Nh, T, Eu\- For V — 2.5 the change in CL is only 14% for more than two orders of magnitude change in the 
column density. For T — 1.9 and 1.5 the change in CL between the two column densities is only a factor of ~ 2.5 and ~ 11 
respectively. The latter still falls short by two orders of magnitude for the change in CL to be commensurate with the change 
in Nh- The second result that is apparent from Fig. 3 is that for T — 2.5 most of the incident energy has been deposited in 
the torus below 10 keV, and the incident X-ray continuum above 10 keV makes an insignificant contribution. On the other 
hand, for T = 1.5, most of the energy is deposited above 10 keV. This is because the shape of the incident spectrum is critical 
on determining whether absorption or Compton scattering dominates the energy losses. Steeper spectra have relatively more 
photons at lower energies than flatter spectra so absorption is correspondingly more important for steeper spectra. In fact, 
Fig. 3 shows that the total integrated energy loss for T = 2.5 is larger for the smaller column density than it is for both 
r = 1.5 and T = 1.9 for the larger column density. In other words, the energy loss for the V = 2.5 Compton-thin case is larger 
than the Compton-thick case for F = 1.5 and F = 1.9. Finally, we also show in Fig. 3 the result (dotted line) for T — 1.9 
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Figure 3. The cumulative fractional energy loss in the torus as a function of energy (i.e. from curves such as those in Fig. 1, integrated 
from 0.05 keV up to an energy, E\ see equation 3). Pairs of curves are shown for a given incident power-law photon index (T), for two 
column densities: 3 X 10 cm" " 2 (red) and 10 25 cm -2 (black). Three pairs of curves are shown, corresponding to T = 1.5,1.9, and 
2.5. The largest value that any of the curves can possibly have is equal to the covering factor, [Af2/(4-7r)] = (a/c), corresponding to all 
the energy intercepted by the torus being captured. These calculations were done for (a/c) = 0.5 and this value is shown by the gray 
dashed line. The dotted curve shows the effect of ignoring Compton scattering for the case of T = 1.9 and Njj = 10 25 cm" 2 (i.e. pure 
absorption). 

when Compton scattering is neglected (i.e. absorption only). It can be seen that for T = 1.9, Compton scattering increases 
the energy loss by ~ 40% compared to the case of absorption only. 

4.1 Dependence on column density 

In this section we examine the explicit dependence of the total integrated fractional energy loss, CL[(a/c), Nh,T, Eu], on Nh- 
Fig. 4 shows this dependence for three values of F (1.5, 1.9, and 2.5) and for each of these three values of V two curves are 
shown, corresponding to two values of the covering factor, (a/c) = 0.1 (dotted) and (a/c) = 1.0 (solid). Note that the quantity 
that is actually plotted in Fig. 4 is CL[(a/c), Nh , T, Eu]/(a/c), the total integrated fractional energy loss per unit covering 
factor. The reason for this is to facilitate a direct comparison of the shapes of the curves for different covering factors. We see 
from Fig. 4 that for T = 2.5, the shape of the total energy loss curves as a function of Nh does not depend on the covering 
factor. For flatter spectra there is some dependence that develops for column densities greater than ~ 10 24 cm -2 but even for 
r = 1.5 the difference between the (a/c) = 0.1 and (a/c) = 1.0 curves is no more than ~ 25% at the highest column density 



The most important result that Fig. 4 shows is, as already discussed above for Fig. 3, the very different behavior of the 
energy loss functions as a function of Nh for different values of F. For the steepest incident X-ray continuum with T = 2.5, 
the sensitivity of CL[(a/c), Nh, I\ Eu] to Nh is less than 20% for a factor of ~ 330 change in Nh- This is because for such a 
steep spectrum most of the energy is in the soft X-ray band, and photons with such low energies are readily absorbed by small 
column densities. Adding more column density then cannot make much difference to the energy loss functions if the bulk of 
the energy has already been deposited in the torus. Fig. 4 shows that even for F = 1.5, the change in CL[(a/c), Nh, T, Eu] 
as JVh changes by a factor of ~ 330 is only just over an order of magnitude. Overall, Fig. 4 shows that the total integrated 
fractional energy loss is actually not very sensitive to Nh for any value of T or (a/c). 

4.2 Dependence on covering factor 

In this section we examine the dependence of the total, integrated fractional energy loss, CL[(a/c), Nh, V, Eu] on the covering 
factor, (a/c). Fig. 5 shows curves of the total integrated fractional energy loss versus (a/c), for five different column densities in 
the range 3 x 10 22 cm -2 to 10 25 cm -2 , calculated for F = 1.9. It can be seen that the relationship between CL[(a/c), Nh,7, Eu] 
and (a/c) for a given value of Nh is very simple. In fact, as might be expected, the relationship is linear for the Nh = 10 24 cm' 2 
curve and for the curves with smaller column densities, but there is little departure from linearity even for the Nh = 10 25 cm -2 
curve. The approximate linearity of CL[(a/c), Nh, T, Eu] as a function of (a/c) even in the Compton-thick regime is true for 
any F in the range relevant to the present study (range 1.5-2.5). The significance of this result is that, since we have already 
seen that the dependence of CL[(a/c), Nh , T, Eu] on iVji can be very weak (much weaker than linear), we can expect the 
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Figure 4. The total 0.05-500 keV integrated energy loss in the torus as a fraction of the total integrated energy in the same bandpass, as 
a function of the torus equatorial column density, Nh- Three pairs of curves are shown, corresponding to photon indices of the incident 
power-law spectrum, T, of 1.5, 1.9, and 2.5. Each pair corresponds to two different values of the covering factor, [Af2/(47r)] = (a/c), of 
0.1 (dotted) and 1.0 (solid). The curves have been divided by the covering factor in order to facilitate a direct comparison of the shape 
of each pair of curves. The functional dependence on Nh does not vary significantly with (a/c). 
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Figure 5. The total 0.05-500 keV integrated energy loss in the torus as a fraction of the total integrated energy in the same bandpass, 
as a function of the covering factor, [Afi/(47r)] = (a/c). Curves are shown for five different column densities, as indicated, for an incident 
power-law photon index of T = 1.9. It can be seen that the total integrated fractional energy loss is approximately a linear function of 
the covering factor. 



covering factor to play a role in the reprocessing of X-rays to infrared emission that is at least, if not more important, than 
that of the column density. 



4.3 Dependence on photon index 

In this section we examine the explicit dependence of the total integrated fractional energy loss, CL[(a/c), Nh, T, Eu], on 
the photon index of the incident power-law continuum, F. Fig. 6 shows this dependence for five different column densities in 
the range 3 x 10 22 cm~ 2 to 10 25 cm~ 2 , for (a/c) = 0.5 (solid curves), as F varies from 1.4 to 2.6. Also shown in Fig. 6 are 
curves for (a/c) = 0.1 and (a/c) — 1.0 for the two highest values of Nh, 5 x 10 24 cm~ 2 (dotted) and 10 25 cm -2 (dashed). The 
energy loss has again been normalized by the covering factor, (a/c), so that the shape of the curves can be directly compared 
for different values of the covering factor. The curves for column densities lower than 5 x 10 24 cm~ 2 in Fig. 6 do not show a 
discernible difference for different values of (a/c). 

We see from Fig. 6 that the total fractional energy loss functions for the smallest column densities are actually the most 
sensitive to V. For the lowest column density (3 x 10 22 cm -2 ), CL[(a/c), Nh , T, Eu] varies by a factor of ~ 18 as T varies 
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Figure 6. The total 0.05-500 kcV integrated energy loss in the torus as a fraction of the total integrated energy in the same bandpass, 
as a function of the photon index of the incident power-law spectrum (T), for five different values of the torus equatorial column density, 
as indicated. The solid lines correspond to calculations with a covering factor, [Af2/(47r)] = (a/c) of 0.5. The dotted lines correspond, for 
Nfj = 5 X 10 24 cm~ 2 , to (a/c) = 0.1 (lower curve), and (a/c) = 1.0 (upper curve). The dashed lines correspond, for Njj = 10 25 cm~ 2 , to 
(a/c) = 0.1 (lower curve), and (a/c) = 1.0 (upper curve). All of the curves have been divided by the covering factor in order to facilitate 
a direct comparison of their shape. 



between 1.4 and 2.6. However, for the highest column density (10 25 cm -2 ), CL[(a/c), Nh, T, Eu] varies by only a factor of 
~ 1.3-1.7 (depending on the covering factor) as T varies between 1.4 and 2.6. The reason for the larger sensitivity for smaller 
column densities is again that the dominant energy loss mechanism is absorption of low-energy photons in that regime. When 
the medium becomes Compton-thick, energy losses due to Compton downscattering dominate and multiple scatterings of 
high-energy photons tend to mitigate the sensitivity to I\ 



5 ENERGY LOSS AS A FRACTION OF THE OBSERVED TO INTRINSIC X-RAY LUMINOSITY 
RATIO 

So far we have seen that of the three critical parameters of the system consisting of the X-ray reprocessor and its illuminating 
continuum (covering factor, JVh, and T), the energy deposited in the medium that is potentially available for reprocessing 
into infrared emission, is least sensitive to Nh • This already does not look promising for supporting the idea that the column 
density of the obscuring matter in AGNs could be constrained by the infrared emission. However, it is not the absolute infrared 
luminosity that is claimed in the literature to be the indicator of column density, but the ratio of the infrared luminosity 
to the observed X-ray luminosity. Since the latter has a strong dependence on the line-of-sight column density, it might be 
supposed that the ratio does in fact provide a good indicator of the column density. In order to address this question, we 
constructed a proxy for the infrared to X-ray luminosity ratio that is due to X-ray energy from the illuminating continuum 
deposited in the obscuring matter. Specifically, using our Monte Carlo results, we calculated the ratio of the total integrated 
fractional energy loss, CL[(a/c), Nh,T, Eu] (see equation 3) to the ratio of the observed to intrinsic X-ray luminosity in the 
2-10 keV band. We call this quantity PX[(a/c), Nh, T], which can be written as 

PX[(a/c),N H ,F] = CL[(a/fN B ,rEv] (4) 

1-^obscrvcd/ -^intrinsic] 2 — 10 kcV 

The quantity CL[(a/c), Nh, T, Eu] represents only the energy potentially available for reprocessing into infrared emission 
so the quantity PX[(a/c), Nh ,F] is really an upper limit for a given set of parameters. However, this does not impact the 
basic test that we want to perform. We need to establish whether the ratio PX[(a/c), Nh, T] is sensitive enough to Nh that 
it could distinguish between a Compton-thin and a Compton-thick AGN, despite the sensitivity of that ratio to V and to 
the covering factor. Deep X-ray survey data on AGNs do not have a sufficiently high signal-to-noise ratio to unambiguously 
determine T and (a/c). Both parameters are highly model-dependent and degenerate with each other and with other model 
parameters for the signal-to-noise ratio that is typical of these deep survey AGN data (e.g., see discussion in Georgantopoulos 
et al. 2009). In the least constrained scenarios, X-ray spectroscopy is not even possible and inferences are made from X-ray 
hardness ratios alone, which carry an even greater degree of degeneracy (e.g., Polletta et al. 2006). We also note that covering 
factors deduced independently from other wavebands (e.g., the dust covering factor) are not necessarily equal to the X-ray 
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N H (10 24 cm" 2 ) 

Figure 7. (a) The ratio of the total 0.05-500 keV integrated fractional energy loss as fraction of the 2—10 keV observed to intrinsic 
luminosity ratio, for a torus viewed edge-on, as a function of the column density. This ratio, PX [(a/c), Nh, T] (see equation 4), is a 
proxy for the ratio of the fraction of the total X-ray luminosity that could be reprocessed into the infrared band, to the 2-10 keV X-ray 
observed to intrinsic luminosity ratio. Curves are shown for two values of the photon index of the incident power-law spectrum, T = 1.9 
(black) and 2.5 (red). For each value of T, the curves are calculated for three values of the covering factor, (a/c) = Af!/(47r) = 0.1 
(dotted), 0.5 (solid), and 1.0 (dashed). The blue crosses joined by the blue horizontal line show that the proxy PX [(a/c), Ng, T] is not 
a good indicator of Nh because, in this example, a Compton-thin source with N H = 2 X 10 23 cm 2 and a Compton-thick source with 
Nh = 1.4 X 10 24 cm -2 have the same value of the proxy. This is because the dependence of the proxy on T and the covering factor can be 
greater than the dependence on Nh . (b) As (a) but this time the 2-10 keV luminosity included a power-law continuum due to scattering 
in an optically-thin zone surrounding the torus, corresponding to 1% of the direct (incident) X-ray continuum. The dependence of the 
proxy PX[(a/c), Nh , T] on Nh is now even weaker. The blue crosses joined by the blue horizontal line show that a Compton-thin source 
with Nh = 2 X 10 23 cm - 2 and a Compton-thick source with Nh = 1.9 X 10 24 cm -2 have the same value of the proxy. 



covering factor. Therefore, in the absence of knowledge of F and the covering factor, if the theoretical ratio PX[(a/c), Nh, T] 
is sufficiently sensitive to distinguish between a Compton-thin and a Compton-thick AGN, then the observed infrared to X-ray 
luminosity ratio might be a useful indicator of the column density. On the other hand, if PX[(a/c), Nh, T] is not sufficiently 
sensitive to distinguish between a Compton-thin and a Compton-thick AGN, then the observed infrared to X-ray luminosity 
ratio is not an indicator of the column density. 

Any additional infrared emission due to dust emission can only serve to reduce the usefulness of the observed infrared 
to X-ray luminosity ratio as an indicator of column density. There is at least one more complication. That is, we know from 
studies of bright AGNs that even in obscured sources the X-ray spectrum often shows an unobscured continuum component 
that is due to optically-thin electron scattering in an extended zone surrounding the central engine. Typically, this component 
can be characterized by a power-law continuum that has a luminosity that is ~ 0.02-5% of the direct (incident) continuum 
(e.g., Turner et al. 1997). However, in X-ray deep survey data the signal-to- noise ratio of the spectra for individual AGN is 
too poor to constrain this component. The resulting spectrum can then appear to have a smaller column density than that 
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which is actually obscuring the primary continuum if the optically-thin scattered continuum is not included in the model (and 
if it is, some assumptions have to usually be made about the pertinent parameters). In calculating the observed to intrinsic 
X-ray luminosity ratio in equation 4 we can include the optically-thin scattered continuum component, with its photon index 
and luminosity relative to the direct (incident) continuum as additional parameters. We note that the observed to intrinsic 
X-ray luminosity ratio depends on the orientation at which the torus is observed. In the following examples we will use an 
edge-on orientation because the the observed to intrinsic X-ray luminosity ratio is most sensitive to A^h when the torus is 
observed edge-on. Since we have already found that the energy losses are not very sensitive to Nh, we are looking for ways that 
PX[(a/c), Nh, r] could recover some sensitivity to Nh and the edge-on orientation is the most appropriate for this purpose. 

In Fig. 7 we show calculations of PX[(a/c), Nh,T] from our Monte Carlo results. The curves in Fig. 7(a) were calculated 
with no optically-thin scattered continuum component and the curves in Fig. 7(b) were calculated with the inclusion of a 
continuum component with the same value of T as the primary continuum but with only 1% of the luminosity of the direct 
(incident) continuum. The curves shown in Fig. 7 were calculated for three values of the covering factor, (a/c) — 0.1 (dotted 
lines), (a/c) = 0.5 (solid lines), and (a/c) = 1.0 (dashed lines), and two values of F (1.9 and 2.5, corresponding to the black 
and red curves respectively). 

We see from Fig. 7 that the sensitivity of PX[(a/c), Nh , T] to the covering factor and V is still very strong and can 
in fact override the weak dependence of PX[(a/c), Nh , T] on Nh- We also see that the effect of including the optically-thin 
scattered continuum, even at the level of 1%, significantly weakens the sensitivity of PX[(a/c), Nh, T] to Nh- For the highest 
column densities, in the Compton-thick regime, we see that PX[(a/c), Nh , T] drops by two orders of magnitude due to the 
optically-thin scattered continuum component, for the same values of V, (a/c), and Nh- 

In both Fig. 7(a) and Fig. 7(b) we illustrate one example in each case of a situation in which a Compton-thin AGN can 
give the same value of the ratio PX[(a/c), Nh, T] as a Compton-thick AGN. The blue crosses joined by the blue horizontal 
lines show that the infrared to X-ray luminosity proxy is not a good indicator of Nh ■ In Fig. 7(a) the example shows that a 
Compton-thin source with Nh — 2 x 10 23 cm~ 2 and a Compton-thick source with Nh = 1.4 x 10 24 cm~ 2 have the same value 
of the proxy. In other words, a Compton-thin AGN with a steep spectrum and a high covering factor could give the same value 
of PX[(a/c), Nh , T] as a Compton-thick AGN with a natter spectrum and a smaller covering factor. In Fig. 7(b) the blue 
crosses joined by the blue horizontal line show that a Compton-thin source with Nh — 2 x 10 23 cm~ 2 and a Compton-thick 
source with Nh = 1.9 x 10 24 cm~ 2 have the same value of the proxy. 

A general conclusion that can be drawn from Fig. 7 is that if we have a large sample of AGNs from a deep X-ray survey 
and examine the distribution of the ratio of observed infrared luminosity to X-ray luminosity, we could say that the AGNs 
with the highest value of that ratio might be Compton-thick candidates and/or heavily contaminated by dust emission. For 
the bulk of the AGNs in the distribution, the infrared luminosity to X-ray luminosity of a given AGN would tell us very 
little about the column density in the absence of information on the intrinsic X-ray continuum, the optically-thin scattered 
continuum, and the covering factor. This realization has already become apparent to Georgantopoulos et al. (2010) who have 
observationally identified some X-ray selected "infrared excess" AGNs that they call DOGs (Dust Obscured Galaxies) to be 
Compton-thin. Our theoretical results explain the reasons why the infrared luminosity to X-ray luminosity in AGNs is not a 
good indicator of the column density. 

Our calculations have assumed a uniform distribution of matter, but the toroidal matter distribution may be clumpy (e.g., 
see Elitzur 2008, and references therein). In the optically-thin limit, a clumpy matter distribution will give the same results for 
the same actual covering factor, except that the covering factor is no longer an indicator of the solid angle subtended by the 
reprocessor at the X-ray source. An additional parameter, a filling factor, would be required to derive a relation between the 
solid angle and the covering factor. In the limit that each clump of matter is Compton thick, the total energy loss to the torus 
would be less than the corresponding uniform matter distribution with the same covering factor and mean column density 
(for the same incident X-ray continuum). This is because more lines-of-sight are available for X-ray photons to escape from 
the surface of a Compton thick clump after one or more scatterings (which may occur in more than one clump) . Therefore the 
overall effect of dumpiness is to make the energy loss even less sensitive to the mean column density of the matter distribution 
because in the optically-thin limit the energy losses are the same but in the Compton-thick limit the losses are less for a 
clumpy distribution compared to a uniform matter distribution. Thus, the principal conclusions of the present paper still hold 
for a clumpy matter distribution. 



6 SUMMARY 

The aim of this paper was to investigate the theoretical foundation of the idea that the observed mid-infrared luminosity 
to X-ray luminosity ratio in AGNs that are found in deep X-ray surveys can be used as a proxy for the column density of 
obscuring matter, which therefore might be useful for identifying Compton-thick candidates. The basis for this is that energy 
lost by the X-ray continuum to the obscuring matter appears as reprocessed infrared emission. In order for the proxy to be 
viable a minimal requirement is that the energy losses due to absorption and Compton scattering in the obscuring matter 
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should be sensitive to the column density and insensitive to other key physical parameters of the system. By means of Monte 
Carlo simulations we found that the energy deposited in the obscuring matter and available for reprocessing into infrared 
emission has a sensitivity to the shape (steepness) of the incident X-ray continuum and the covering factor of the intercepted 
material that is far greater than the sensitivity to the column density. As a result, we found that the observed infrared to 
X-ray luminosity ratio for a Compton-thin AGN could be just as large as that for a Compton-thick AGN. The signal-to-noise 
ratio of the X-ray spectra of AGNs currently found in deep X-ray surveys is not sufficiently high to constrain the intrinsic 
X-ray spectrum and the covering factor of the reprocessor well enough to enable the infrared to X-ray luminosity ratio to be 
used as a reliable indicator of column density. If the infrared emission has a significant contribution from dust emission, for 
example due to starburst activity, the reliability of that indicator can only be worse. Our conclusion regarding the difficulty 
in the use of the mid-infrared to X-ray luminosity ratio for distinguishing between Compton-thin and Compton-thick AGNs 
is independent of the geometry of the obscuring matter in AGNs. It is also independent of the angular distribution of the 
intrinsic X-ray continuum emission. This is because the dominant energy loss mechanism changes from absorption to Compton 
scattering if the slope of the X-ray continuum is changed from having a power-law photon index of 2.5 to 1.5 (the typical 
range in AGNs), and this fact is independent of geometry. Further, the energy deposited in the obscuring matter trivially 
depends on the covering factor, and this fact is also independent of geometry. 
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